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Abstract 

We study the stability and breakup of liquid bridges with a free 
contact line on a surface with contact-angle hysteresis under zero- 
gravity conditions. Theoretical predictions of the stability limits are 
validated by experimental measurements. Experiments are conducted 
in a water-methanol-silicon oil system where the gravity force is off¬ 
set by buoyancy. We highlight cases where stability is lost during 
the transition from a pinned-pinned to pinned-free interface when the 
receding contact angle is approached—rather than a critical state, in¬ 
dicating that the breakup length is not always associated with the 
static maximum-length stability limit. We demonstrate that the dy¬ 
namic contact angle controls the contact-line radius following stability 
loss, and that interface evolution following stability loss can increase 
the dispensed-drop size if the contact angle is fixed. 


1 Introduction 


Studying the stability and dynamics of liquid bridges is motivated by a broad 


range of applications, including crystal growth in microgravity (Garcia Ve¬ 


larde, 1988), surface patterning, nano-printing, and nano-lithography (Salaita 
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Quantifying liquid-bridge and jet breakup upon stability loss dates to the 
works of Plateau (1873) and Rayleigh ( 1879ap> ). While early investigations 
on crystal growth and purification focused on determining the minimum liq¬ 


uid volume that can be held between circular discs (Martinez & Perales 


1986 Meseguer et al., 1995), the drop-size distribution following breakup 


is of prime interest in contact-drop dispensing and liquid-transfer applica¬ 


tions (Dodds et al. 2009 Lutfurakhmanov et al. 2010 Chen et al . 2014) 


Minimizing the dispensed-drop size relative to the needle diameter is central 


to surface patterning based on direct-write lithographic techniques (Piner 


et al ., 

1999 

Huo et al., 

2008) 


Recent studies on contact-drop dispensing have shown that the deposited 
drop size is influenced by the needle retraction speed, needle-tip size, sur¬ 


face characteristics, and dispensing control parameters (Qian et al., 2009 


Qian & Breuer, 2011). Interestingly, the deposited drop volume in pressure- 


controlled and volume-controlled dispensing behave differently with the nee¬ 
dle retraction speed. Faster retraction reduces the drop size to a minimum 
and monotonically increases the drop size in pressure-controlled and volume- 
controlled dispensing, respectively. Three regimes were experimentally iden¬ 
tified with respect to the retraction speed U n for the pressure-controlled case. 
In the first two, U n <C u w , where u w is the capillary-wave speed; the con¬ 
tact line is advancing in the first and stationary in the second, and the 
drop size scales as f/T 1 ^ 2 ; the third corresponds to fast retraction speeds 
( U n /u w ~ O(10~ 2 )) where the dynamics dramatically change, and the drop 
size does not scale with U n as a simple power law. Here, the drop size is 
almost two orders of magnitude smaller than in the first two regimes, which 


Qian et al. (2009) attributed to a fast receding contact line with a speed ap¬ 


proaching u w . However, in volume-controlled deposition, the dispensed-drop 
size did not exhibit the same sensitivity to the needle retraction speed in 
the parameter range studied by Qian & Breuer (2011). Thus, our study is 
motivated, in part, by seeking to answer whether it is possible to influence— 
by purely geometric means—the dynamics in volume-controlled deposition, 
so that one may achieve comparable sensitivity as in the pressure-controlled 
case, to achieve small-drop deposition. 
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Liquid transfer in contact-drop dispensing is generally accomplished by 
the breakup of a drop bridging a sharp-edged needle and a flat substrate. 
The bridge is pinned to the needle edge, and may have a free or pinned 
contact line with the substrate, depending on the surface wettability. To 
help distinguish the effects of geometric parameters on the dispensed-drop 
size from dynamics ones, Akbari et al. (2015c) studied the static stability of 
liquid bridges with free-pinned contact lines, showing that the contact-line 
radius and bridge-neck position—key parameters determining the dispensed- 
drop size—can be controlled by purely geometric means. Constructing the 
stability region at fixed contact angles, they also showed that the free contact 
line has symmetry-breaking and destabilizing effects (Akbari et al ., 2015o[c). 

Experimental studies on the statics of liquid bridges between equal cir¬ 
cular discs are extensive. Using neutral-buoyancy experiments, Sanz fc Mar¬ 


tinez (1983) ascertained the minimum-volume stability limit in the slender¬ 


ness range 0 < A < 6. Russo & Steen (1986) determined the maximum- 


volume stability limit in a similar set-up, showing that axisymmetric liquid 
bridges non-axisymmetrically bulge when their interface is tangent to the 
discs. The experiments of Slobozhanin et al. (1997) provided further in¬ 
sights on this stability limit, ffere, the stability limit corresponds to pitch- 
fork bifurcations. Moreover, they showed that, above (below) the slen¬ 
derness A ~ 0.4946, liquid bridges continuously (abruptly) bulge into a 
non-axisymmetric shape. Other studies considered the effect of gravity on 


the stability limits of axisymmetric ( 

Bezdenejnykh et al. 

1992 

) and non- 

axisymmetric liquid bridges ( 

Bezdenejnykh et al., 

1999 

) between equal discs. 


Emphasizing the destabilizing effect of gravity on nearly cylindrical liquid 
bridges, Lowry & Steen (1994) experimentally demonstrated that subjecting 
liquid bridges to an external laminar flow suppresses interfacial disturbances, 
thereby stretching the stability limit beyond that of static bridges. 

Surface imperfections (e.g., heterogeneity and roughness) complicate the 
equilibrium of gas-liquid-solid contact lines on real surfaces compared to ideal 


Joanny & de Gennes 

1984 

Qian & Breuer 

2011 

Chen et al. 


2013). Contact lines remain pinned on real surfaces as long as the equilibrium 


contact angle is between the receding and advancing contact angles, and 
otherwise freely move (Gao & McCarthy, 2006). The receding contact angle 
is a key geometric parameter that affects the dispensed-drop size when the 


contact line is free (Akbari et al.. 2015c Qian & Breuer, 2011). Chen et al 


(2013) experimentally and numerically studied the effect of contact-angle 
hysteresis on the the evolution and adhesion force of liquid bridges with 
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two free contact lines. Similarly to Qian & Breuer (2011), experiments were 
conducted in a liquid-gas system where the gravity effect is alleviated by small 
bridge dimensions. Using similar experiments, Chen et al. (2014) examined 


the bridge breakup, showing that the liquid transfer ratio is correlated with 
the difference between the receding contact angles on the plates. 

When stretching liquid bridges, there is a critical length, identified with 
the static stability limit (Qian & Breuer, 2011), where the dynamics dramat¬ 
ically change. Below this length, the dynamics are slow and the interface 
evolves quasi-statically. However, the dynamics are fast above the critical 
length where fluid viscosity and inertia become significant as the interface 
nears pinch-off, and the liquid and interface velocities approach u w . Because 
the dynamics in the pinch-off phase are much faster than in the quasi-static 
phase (Eggers, 1997), the pinch-off length is commonly approximated as the 
static maximum-heigh stability limit in the literature (Qian & Breuer 2011[ 


Chen et al. 2014) 


Studying the relationship between the contact-line constraints and the 
stability limits of liquid bridges, Akbari et al. (2015c) demonstrated that 
pinned-pinned liquid bridges can be stretch beyond the maximum-length 
stability limit of pinned-free bridges for the same liquid volume. Conse¬ 
quently, unstable states of the latter can be accessed during drop deposition 
if surface imperfections constrain the disturbances at the contact line. This 
has significant implications for the breakup dynamics and liquid transfer 
on (real) surfaces with contact-angle hysteresis. Depending on the surface 
wettability and drop volume, the contact angle may fall between or outside 
the receding and advancing contact angles when stretching liquid bridges. 
Hence, the bridge can undergo transitions from pinned-pinned to pinned-free 
contact lines (and vice versa) during stretching. These complications raise 
new, non-trivial questions as to how the stability limits and dispensed-drop 
volume are influenced by contact-angle hysteresis, and which stability limit 
(with respect to pinned-pinned or pinned-free disturbances) determines the 
breakup length. 

In this paper, we address the forgoing questions by studying pinned- 
pinned to pinned-free transitions and their respective stability limits during 
drop deposition on surfaces with contact-angle hysteresis. In particular, we 
show that, contrary to the common notion, there are cases where liquid 
bridges do not break at a critical state (he., the pinned-pinned or pinned- 
free stability limit). This observation has not been reported in the literature, 
as far as we are aware. Moreover, we experimentally verify the theoretical 
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predictions of the maximum- and minimum-slenderness stability limits (Ak- 


bari et al. 2015c) when the contact line is free at breakup. The stability 


limits of liquid bridges with pinned and moving contact lines are compared, 
demonstrating the destabilizing effect of free contact lines. To simulate zero 
gravity using the density matching technique, liquid bridges of silicon oil 
were formed in a water-methanol solution. The contact-angle effect was then 
studied by adding surfactant to the aqueous phase. 


2 Materials and methods 

Experiments were performed in a cubic Plateau tank under neutrally buoyant 
conditions (Fig. [Tj) . Silicon oil (5 cSt, Sigma Aldrich) with specific gravity 
0.92 was used as the bridge in a water-methanol solution (volumetric mix¬ 
ing ratio 42:58) bath. The composition of the bath solution was adjusted 
so that its density matched that of silicon oil at the experiment tempera¬ 
ture (~ 20°C). Using a microsyringe, a drop with a prescribed volume in 
the range 5-50 /A was deposited onto a plastic coverslip (Fischer Scientific), 
which had been soaked in a 0.1 M hydrochloric acid solution, rinsed with DI 
water, and placed in the tank. A bridge was produced by gently pressing a 
needle with tip diameter 1.5 mm into the drop. The needle was mounted on 
a one-dimensional vertical translation stage to control the bridge length, and 
the tank was placed on a two-dimensional positioning stage to align the drop 
and needle centers before contact, ensuring that the bridge is axisymmet- 
ric. The maximum (minimum) slenderness stability limit was ascertained by 
stretching (squeezing) the bridge until the bridge ruptured (bulged asymmet¬ 
rically). A CCD camera (Prosilica GX1050, Allied Vision) with a 5x lens 
(Nikon GMicro-NIKKOR) was used to record the bridge dynamics. Images 
were analyzed using an in-house Matlab script. The bridge contact angle 
with the coverslip was adjusted by changing the interfacial tensions in the 
system using the anionic surfactant sodium dodecyl sulfate (SDS) (Sigma 
Aldrich) at concentrations in the range 0-10 g 1 _1 . 

3 General behaviour of liquid bridges 

Since the is was hollow with a sharp edged tip, bridges were always pinned 
to the needle. After a drop was deposited onto a coverslip, it was squeezed 
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Figure 1: Schematic of the experimental setup. 


in 0.01 inch steps to reach the minimum-slenderness stability limit, at which 
it bulged asymmetrically; bridges were imaged at each step. The maximum- 
slenderness stability limit was similarly measured by stretching bridges until 
rupture. Figure [2] shows a typical sequence during squeezing and stretching 
of a 20 pi drop. In the rotund limit, non-axisymmetric drop deformations 
displace the contact line, leading to different results upon repeating the ex¬ 
periment. In the slender limit, the bridge breaks into two primary drops, 
leaving several satellite drops suspended in the bath. Without SDS, the con¬ 
tact line moved only at breakup for small drop volumes (less than 10 pi), 
and was otherwise pinned at larger volumes. Since the emphasis in this work 
is on the role of moving contact lines, SDS was added to the bath to reduce 
the receding contact angle (see Fig. [3]) of the bridge on the coverslip. 

4 Theory 

Consider a liquid of volume v bridging a circular disk with radius Rq and a 
large plate. The disc and plate are separated by a distance h , as shown in 
Fig. [3j The bridge is pinned to the disc and is free to slide horizontally on 
the plate. A small Bond number (Bo -C 1) is achieved by density matching, 
so the gravity force is negligible. Consequently, there is a constant pressure 
differential between the non-hydrostatic pressure of the bridge pi and the 
surrounding fluid p g . Here, the surface tension between the phases i and j is 
denoted yjj with F^- the corresponding interfacial surface area. The contact 
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Figure 2: Representative stretching (bottom) and squeezing (top) sequences 
(20 //I drop). 


and dihedral angles that the interface T g i forms with the plate and disc are 
denoted 9 C and 9 d , respectively. The cylindrical volume V = v/('kRqK), 
scaled volume v* = v/(4irRl/3), scaled pressure (mean curvature) Q = qRo, 
and slenderness A = h/Ro are the dimensionless parameters with which the 
liquid bridges are specified. Note that q = (p g — Pi)/^ g i measures the non¬ 
hydrostatic pressure differential (Myshkis et al, 1987). Solving the Young- 


Laplace equation furnishes the equilibrium meridian curve 


p(r) = y/l + a 2 + 2a cost 

= So t 


(i) 


where a = p(0) — 1 and r is the mean-curvature-scaled arclength (Akbari 


et al, 2015c). The scaled arclength at the hole edge £ and the contact line 


£ are denoted r 0 and ri, respectively. Akbari et al. (2015c) showed that 


critical surfaces at the minimum-slenderness stability limit, corresponding to 
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Figure 3: Weightless liquid bridge; (a) schematic and (b) coordinate system 
with meridian curve parametrization. 


the upper boundary of the stability region, are nodoids with 9 d = 0, which 
are well approximated by 

V = 1 + - sec 4 (6 ) c /2)(7r — 6 C + cos6 ) c sin6 l c )A 

- sec 8 (6 l c /2)[—97 + 24 (t r - 9 C ) 2 - 136 cos 9 C - 32cos(2 9 C ) ^ 

384 

+ 8 cos^c) + cos(46h) + 24 ( 7 t — 9 C ) sin(26 ) c )] A 2 + 0(A 3 ) 


in V versus A stability diagrams. 

Note that pinned-free and pinned-pinned liquid bridges are, respectively, 
specified by p = (A, V, 9 C ) and (A, V, K ), where K is the ratio of the lower to 
upper contact-line radii. When stretching the bridge, K varies at fixed 9 C in 
pinned-free bridges, and 9 C varies at fixed K in pinned-pinned bridges. Fur¬ 
thermore, we determine stability along equilibrium branches using Myshkis’s 
variational method, detailed by Akbari et al. (2015c). 


5 Feature extraction 


The bridge interface with the bath solution was extracted using a gradient- 


based edge detection method with a Gaussian optimal smoothing filter (Marr 


& Hildreth. 1980). In this method, pixels on an interface are identified by 


finding maxima in the first directional derivative of intensity; or, equiva¬ 
lently, seeking zero-crossings in the second directional derivative. Derivatives 
were taken along normals to interfaces using high-order (8-10 points) central 

























Figure 4: Fitting theoretical meridian curves to bridges boundaries in stretch¬ 
ing (bottom) and squeezing (top) sequences (20 /il drop). 


schemes. Then, the analytical solution of the bridge meridian curve, given by 
Eq. Q, was fitted to the extracted interfaces. Here, the unknown parameters 
(Q c , To, Ti, a) were determined by minimizing the root-mean-squared normal 
distances between the extracted interface pixels and the theoretical merid¬ 
ian curve. Figure [4] shows typical results of the image-processing script in 
stretching and squeezing experiments. 


6 Results and discussion 


6.1 Surfactant effect 


As previously stated, for large drops, the contact angle 9 C remains smaller 
than the receding contact angle 6 r during stretching. Thus, the contact 
line is pinned to the coverslip at breakup. To assess the stability limits 
of liquid bridges with a free contact line over a wider range of drop vol¬ 
umes, the contact angle was reduced by adding SDS to the bath solution. 
The critical micelle concentration (CMC) of SDS in pure water at 25°C is 
~ 2.36 g F 1 (Mukerjee & Mysels, 1971). Previous measurements in the lit¬ 


erature have shown that the air-water surface tension exhibits no minimum 
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Figure 5: The surfactant-concentration effect on the sessile-drop contact an¬ 
gle 9 s d at drop volumes 5 /if (o), 10 /A (A), 15 [A (□), 20 /A (x). Dashed line 
indicates the critical micelle concentration of the surfactant in pure water at 
25°C (iMukerjee & Myselsl Il97lb. 


near the CMC ( 

Lucassen-Reynders et al. 

1981). According to Young’s equa- 

tion ( 

Myshkis et al, 

1987 

), this implies that the contact angle also does not 


exhibit a minimum if the air-solid and water-solid surface tensions remain 
constant. Furthermore, previous reports in the literature indicate that, un¬ 
like the advancing contact angle, the receding contact angle can be sensitive 
to the drop volume (Drelich et al, 1996). Therefore, we examine how the 


contact angle varies with the SDS concentration at various drop volumes, in 
the silicon oil-water-methanol system, to determine the surfactant concen¬ 
tration at which the contact angle is minimum for all volumes. The contact 


angle was measured using the sessile-drop method (Bachmann et al, 2000). 


Figure [5] shows the surfactant effect on the contact angle. The contact 
angle decreases almost linearly around the CMC and below ~ 6 g F 1 for all 
drop volumes. At higher concentrations, the relationship is nonlinear. Never¬ 
theless, the smallest value of the contact angle in the range 2-10 g D 1 occurs 
at 10 g l -1 for all volumes, except 5 /A. Moreover, the contact angle for larger 
drops is affected more by the surfactant at this concentration. Since larger 
drops tend to have a pinned contact line at breakup more often than smaller 
drops, only the stability-limit results for experiments where the SDS concen¬ 
tration is 10 g l -1 are reported for all drop volumes. At this concentration, 
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Figure 6: Comparison of the theoretical prediction and experimental mea¬ 
surement of the stability limits with drop volume v ~ 5 /if, receding contact 
angle 9 r ~ 110°, advancing contact angle 9 a « 70°, and without surfactant 
(right). An image sequence of the bridge evolution corresponding to the data 
points (left). Dashed blue and black lines respectively indicate the constant-u 
isocontour at the dispensed drop volume and the maximum-volume stability 
limit estimated by Eq. (J2j) at 9 C — 70°. Labels denote the contact angle in 
degrees. 


the bridge contact line with the coverslip moved at all drop volumes. 

6.2 Stability limits 

The maximum-slenderness stability limit was determined by stepwise quasi¬ 
static stretching of a liquid bridge with fixed volume. Similarly, the minimum- 
slenderness stability limit was determined by stepwise quasi-static squeezing 
of a bridge. In all experiments, the contact line was either pinned or receding 
when stretching, and always advancing when squeezing. Therefore, when the 
contact line is moving, the receding contact angle 9 r is the relevant contact 
angle in stretching, and the advancing contact angle 9 a is the relevant con¬ 
tact angle in squeezing. Figure [6] shows an image sequence during stretching 
and squeezing of a 5 (A drop. Here, no surfactant was added to the bath, 
and the contact line moved during stretching and squeezing. The receding 
contact angle at breakup was measured 9 r ~ 110°; thus, the corresponding 
data point (far right) in the stability diagram is expected to fall between the 
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2015c) for 9 C = 90° 

as demonstrated in Fig. |6j The advancing contact angle was mea- 


lower boundary of the stability region (Akbari et al. 
and 120 

sured 9 a ~ 70°; here, the minimum-slenderness stability limit is estimated 
by Eq. (|2| and then compared with the measured value. As shown in Fig. [6j 
experimental measurements of the stability limits are in good agreement with 
the theoretical predictions of Akbari et al. (2015c). Note that, because the 
needle is hollow, part of the initial drop volume is driven into the needle in 
squeezing experiments, so the bridge volume does not reflect the initial drop 
volume; moreover, data points deviate more from the constant-r? isocontour 
corresponding to the initial dispensed volume (dashed blue line) near the 
upper boundary of the stability region. 

Stretching and squeezing experiments were conducted in the range v = 5- 
20 //I with SDS added to the bath. At 10 g l -1 SDS, the advancing and 
receding contact angles drop to 9 a ~ 0-5° and 9 r ~ 75-95°. Reasonable 
agreement is observed between the experimental measurements of the stabil¬ 
ity limits and theoretical predictions of Akbari et al. (2015c) (see Fig. [7]). 

Figure [8] shows the effect of a free contact line on the maximum-slenderness 
stability limit. Here, a stretching experiment was conducted on two cover- 
slips using a 20 /A drop with 10 g 1 _1 SDS. These coverslips exhibited slightly 
different contact angles at breakup, presumably due to different surface char¬ 
acteristics, so that 9 C was below 9 r at breakup on one (Fig. [8j left panel) and 
9 C reached 9 r before breakup on the other (Fig. [8j right panel); consequently, 
the contact line was pinned on the former and free on the latter, and the 
slendernesses at breakup were measured A^ ps 5.99 and A& ps 4.94, respec¬ 
tively. This pa 20% decrease in the breakup length reflects the destabilizing 
effect of a free contact line, as theoretically predicted for static catenoidal 
and cylindrical liquid bridges (Akbari et al. 2015a). 

Figure [9] shows the radius of the meniscus contact line R\ and contact 
angle 9 C in the experiments depicted in Fig. [HJ For the liquid bridge shown 
in the left panel of Fig. [8j the contact line is pinned, and the contact angle 
varies during stretching. Here, the contact angle remains below the receding 
contact angle during the entire experiment. By contrast, for the liquid bridge 
shown in the right panel, the contact angle reaches the receding contact angle 
before breakup, and the contact line recedes during stretching. Here, the 
contact-line motion when stretching the bridge is significantly larger than in 
the left panel, indicating that the contact line on the coverslip shown in the 
left panel of Fig. [8] is more constrained than that in the right panel. 

Figure [9] also illustrates a qualitative difference between systems with and 
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without contact-angle hysteresis. Advancing and receding experiments were 
performed on separate coverslips for the bridge with pinned contact lines 
(triangles). Therefore, the bridge does not experience hysteresis, and the 
contact line and contact angle vary continuously from stretching to squeezing. 
By contrast, squeezing and stretching were consecutively performed on the 
same coverslip for the bridge with a free contact line (circles). Here, the 
contact line retreats on a surface that is covered by the silicon oil during 
squeezing, which affects q / sg as the bridge is stretched. Consequently, the 
receding contact angle is different at a given contact-line position during 
stretching and squeezing, so that the contact line radius follows a different 
path when the bridge is stretched to its initial length and beyond. 


6.3 Contact-angle hysteresis effect 


In this section, we focus on data for the bridge with a free contact line at 
breakup (circles) in Fig. [9] and examine how the breakup length is associ¬ 
ated with the maximum-length stability limit with respect to pinned-pinned 


and pinned-free disturbances. Figure 10 shows the bridge-evolution images, 
comparing the theoretical prediction of the contact-line radius and measured 
values. Reasonable agreement is observed between the measurements and 
theoretical predictions. The last image H was recorded between stability 
loss and breakup. Here, the coverslip exhibited an advancing contact angle 
9 a ~ 3° and two distinct receding contact angles 6 r \ ~ 30° and 0 r2 ~ 75°. 
As previously stated, this can be attributed to changes in the coverslip in¬ 
terfacial tension with the bath solution upon retreating the contact line on 
a surface that is already covered by silicon oil. 

The contact line was free during the entire squeezing experiment, and the 
contact angle remained almost fixed at 6 C = 9 a . However, during stretching, 
the contact line was pinned (K & 0.7699) when 6 a < 9 C < 9 r i, free while 
the contact angle was almost fixed at 9 C = 9 r i, and pinned (K ~ 0.5271) 
when 9 rl < 6 C < 0 r2 . The trajectory P\P 2 indicates the equilibrium solution 
corresponding to the squeezing, whereas P 2 P 3 , P 3 P 4 , and P 4 P 5 correspond to 
the foregoing pinned-pinned, pinned-free, and pinned-pinned phases of the 
stretching experiment. The point P 5 corresponds to a state where 9 C = 9 r2 , 
indicating a transition from pinned-pinned bridges at K = 0.5271 to pinned- 
free bridges at 9 C = 75°. Pinned-free bridges along P 3 P 4 are not stretched 
beyond their stability limit (thin dashed-dotted line), niether are pinned- 
pinned bridges along P 4 P 5 beyond theirs (thin dashed line). Note that the 
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slenderness at P 5 is also well below the pinned-free stability limit (thick 
dashed line) where the foregoing transition occurs. However, experimental 
data indicate that the bridge loses stability during the transition, and it does 
not correspond to the pinned-pinned or pinned-free stability limit. 

To understand why stability is not lost at a critical state, we locate the 
stretching and squeezing trajectories of Fig. [10] on equilibrium branches in 


Figs. 11 and 12 Here, the bridge follows stable equilibrium branches along 
the entire P 1 P 5 trajectory and before the transition at P 5 . We note that 
the bridge lies on a stable branch with respect to pinned-pinned disturbance 
(see Fig. [I2|j) before the turing point. Upon transition at 9 C = 6 b 2 , the 
bridge is exposed to pinned-free disturbances (P 5 in Fig. 11 3 ). Although the 


slenderness is smaller than the maximum-slenderness stability limit (at the 
turning point), this state lies on an unstable branch, where the bridge loses 
stability to pinned-free perturbations. This implies that, it is possible for 
liquid bridges to break during the pinned-pinned to pinned-free transition 
rather than a critical state at a turning point in systems with contact-angle 
hysteresis; therefore, the breakup length is not always associated with the 
static maximum-length stability limit. 

These observations also have significant implications for the breakup dy¬ 
namics and dispensed-drop volume. Recall, pinned-pinned liquid bridges are 
more stable than pinned-free ones. Consequently, the foregoing transition 
provides access to unstable states of pinned-free bridges, far from critical 
sates. Before the transition, the contact line is pinned, and, depending on 
the bridge volume and receding contact angle, can be stretched to a state 
that is highly unstable to pinned-free perturbations. Once the receding con¬ 
tact angle is reached, the constraint at the contact line is relaxed, and the 
interface is exposed to a larger set of perturbations. This leads to a dramatic 
stability loss at a point that does not coincide with a critical state (P5 in 
Figs. 12 d and 11:). Note that the energy barrier (potential well) disappears 
at critical states, and the instability margirQ generally grows as an unstable 
state moves farther away from its critical state along the respective equilib- 


lr The notion of the instability margin in this paper is the stability-margin counterpart 
of Slobozhanin et al. (2002) for unstable bridges. Note that liquid bridges break into 


several primary and satellite drops upon stability loss at the maximum-length stability 


limit (Meseguer et al. 1995). Thus, the chain of drops arising upon breakup is the most 


stable (having the deepest potential well) state that is dynamically accessible to unstable 
bridges. Accordingly, the instability margin is defined as the potential-energy difference 
between the most stable and unstable states. 
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rium branch (Myshkis, 1965 Slobozhanin et al. 2002). Potential energy of 
the instability margin can be transformed to kinetic energy upon stability 
loss, acting as the deriving force for the near-singularity dynamics (EggersJ 


1997). 


An appreciable difference between the breakup dynamics upon stabil¬ 
ity loss at a critical state and at an unstable state away from its critical 
state is expected. At a critical state (e.g., turning point), the energy bar¬ 
rier disappears in the direction of the critical perturbation, so the interface 
accelerates in the same direction with a kinetic energy that is proportional 
to the disturbance magnitude. However, at an unstable state away from its 
critical state, the instability margin significantly amplifies the critical per¬ 
turbation, leading to more dramatic dynamics and a significant impact on 
the dispensed-drop size. As previously stated, Qian et al. (2009) reported 


small drop sizes in a pressure-controlled deposition due to fast dynamics near 
the contact line. Here, achieving a fast-receding contact line is assisted by 
the withdrawal of the liquid near the plate, which is less significant in the 
volume-controlled case (Qian & Breuer, 2011). However, breakup at an un¬ 
stable state with a large instability margin can greatly influence the dynam¬ 
ics in volume-controlled deposition, potentially favouring smaller deposited 
drops. Further studies are required to examine the possibility of stability loss 
at unstable states with large instability margins, and their ensuing dynamics. 

We further elaborate on the behaviour of the contact-line radius after 
stability loss by constructing the equilibrium branches at contact angles near 
0 r 2 . The scaled form of the slender-jet approximation (Eggers & Dupont, 


1994) suggests that, except very close to the singularity where the capillary 
number scales as Ca ~ 0(1), the bridge prohle can be reasonably approxi¬ 
mated by the Young-Laplace equation. This approximation is expected—also 
experimentally shown by Qian & Breuer (2011)—to be accurate for volume- 
controlled stretching since the liquid velocity inside the bridge is restricted 
by the volume constraint. Therefore, the bridge must evolve along unstable 
equilibrium branches at fixed volume. We apply this approximation as a 
guide to investigate the bridge evolution after stability loss and away from 
the pinch-off (e.g., bridge H in Fig. 10) and provide a better understanding 
of the relationship between the receding contact angle and dispensed-drop 
volume. 

Akbari et al. (2015c) showed that, at fixed 6 C , liquid bridges with a free 


contact line exhibit a transcritical bifurcation at a point along the lower 
boundary of the stability region. Transcritical bifurcations were represented 
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by equilibrium branches at fixed 9 C and A in pressure versus volume diagrams. 
Figure [13] shows equilibrium branches in the vicinity of a transcritical bifurca¬ 
tion in R\/Rq versus A diagram at fixed v* and 9 C . As discussed for Figs. 12 d 


and 115. the bridge after the pinned-pinned to pinned free transition at P 5 
lies on an unstable segment of the primary branch at 9 C = 75°. Upon sta¬ 
bility loss, the dynamic contact angle increases from « 75 to 81° with the 
bridge moving from the primary branch of 9 C = 75° at P 5 to the secondary 
branch of 9 C = 81° at H. Here, the opposite behaviour of Pi with A along 
stable and unstable branches is notable: the contact-line radius decreases (in¬ 
creases) during stretching along stable (unstable) branches, favouring small 
(large) drops. This implies that evolution upon the loss of stability at a 
fixed-contact-angle leads to larger drops; thus, to achieve smaller drops, the 
dynamic contact angle must increase following stability loss, upon which the 
contact line accelerates toward pinch-off. Therefore, there are cases where a 
dynamic contact angle that deviates from the receding contact angle plays a 
key role in determining whether the contact line expands or contracts. 


7 Concluding remarks 

We have experimentally and theoretically studied the stability and breakup 
of weightless liquid bridges on surfaces with contact-angle hysteresis. Ex¬ 
periments were performed in a Plateau tank where the effect of gravity was 
alleviated by density matching. To achieve free and pinned contact lines, the 
contact angle was adjusted by adding SDS to the bath. For liquid bridges 
with a free contact line, experimental measurements validated the theoret¬ 
ical predictions of the stability limits. At fixed volume, liquid bridges with 
a free contact line exhibit a larger breakup length than those with a pinned 
contact line, demonstrating the destabilizing effect of a free contact line, as 


theoretically predicted by Akbari et al. (2015c). 


We examined the effect of contact-angle hysteresis on the maximum- 
length stability limit, showing that the breakup length can not always be 
associated with the static stability limit. Depending on the drop volume and 
receding contact angle, liquid bridges may lose stability during the pinned- 
pinned to pinned-free transition at an unstable state away from its critical 
state. This has significant implications for the dynamics following stability 
loss. Unstable states that are far from critical states generally have large 
instability margins, which can transform potential energy to kinetic energy 
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upon stability loss, having a significant impact on the breakup dynamics and 
dispensed-drop size. Constructing equilibrium branches in the vicinity of the 
receding contact angle revealed that a complex interplay between the dy¬ 
namic contact angle (determined by its speed), receding contact angle when 
losing stability, and bridge volume determine whether the contact line is ex¬ 
panding or contracting upon stability loss. Furthermore, we showed (by one 
example) that the contact-line motion upon stability loss at fixed contact 
angle can be advancing, thus disfavouring small-drop deposition. 

More comprehensive computational investigations of the stability and 
breakup of liquid bridges will hopefully provide deeper insights into the re¬ 
lationship between the contact-angle hysteresis and dispensed-drop size, and 
on the effectiveness of surface hydrophobization (to modulate the contact 
angle) for micro-deposition. 
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(a) ( b ) 




Figure 7: Same as Fig. [6j but with SDS concentration 10 g 1 1 at drop 
volumes (a) v = 10 //1, (b) v = 12.5 /il, (c) v = 15 //l, and (d) u = 17.5 //l. 
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Figure 8: Contact-line effect on the breakup length of liquid bridges. Stretch¬ 
ing a 20 /j 1 drop with 10 g l -1 SDS in the bath, producing a pinned contact 
line with 9 C & 84° (left) and a free contact line with 9 C « 81° (right) at 
breakup. 
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Figure 9: Squeezing (filled markers) and stretching (open markers) of a 20 fi 1 
drop with 10 g l -1 SDS, corresponding to the experiments shown in Fig. [ 8 j 
The radius of the meniscus contact line Ri (left) and contact angle 9 C (right) 
versus slenderness are plotted when the bridge contact-line on the coverslip 
at breakup is free (o, right panel in Fig. [ 8 ]) and pinned (A, left panel in 
Fig. Q . 
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Figure 10: Comparison of experimental measurements (open circles) and the¬ 
oretical predictions (solid lines) of the contact-line radius during the stretch¬ 
ing and squeezing of a 20 /d drop (v* ~ 7.171) on a substrate with an ad¬ 
vancing {9 a ~ 3°) and two receding (9 r i ~ 30°, 9 r 2 ~ 75°) contact angle(s). 
Vertical lines indicate the minimum-high stability limit at 9 C = 3° (thick 
dashed-dotted) and maximum-high stability limits at 9 C = 30° (thin dashed- 
dotted), 9 C = 75° (thick dashed), and K = 0.5271 (thin dashed), where the 
contact line is free for the first three and pinned for the last. Equilibrium 
states are computed at fixed contact angle along PiP 2 ( 9 C = 5°) and P3P4 
( 9 C = 30°), and at fixed contact-line radius along P2P3 [K = 0.7699) and 
P4P5 (K = 0.5271). 
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Figure 11: Free-contact line equilibrium branches of a hxed-volume (v* = 
7.171) liquid bridge, indicating stable (solid) and unstable (dashed) states at 
(a) 9 C = 3°, (b) 9 C = 30°, and (c) 9 C = 75°. The terminal points Pi_ 5 are 


shown to identify the stability of equilibrium trajectories in Fig. 10 
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Figure 12: Pinned-contact line equilibrium branches of a fixed-volume (v* = 
7.171) liquid bridge, indicating stable (solid) and unstable (dashed) states at 
(a) K = 0.7699 and (b) K = 0.5271. The terminal points P 2 _5 are shown to 


identify the stability of equilibrium trajectories in Fig. 10 



Figure 13: Free-contact line equilibrium branches of a fixed-volume (v* = 
7.171) liquid bridge in the vicinity of the transcritical bifurcation, indicating 
stable (solid) and unstable (dashed) states. Numeric labels denote the con¬ 
tact angle 9 C in degrees. Circles indicate the states at G, H, and P 5 as the 


bridge approaches breakup in Fig. 10 
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